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Abstract: We introduce field theory techniques through which the deconfinement transition 
of four-dimensional Yang-Mills theory can be moved to a semi-classical domain where it 
becomes calculable using two-dimensional field theory. We achieve this through a double- 
trace deformation of toroidally compactified Yang-Mills theory on R 2 xS[xSj. At large N, 
fixed-L, and arbitrary /3, the thermodynamics of the deformed theory is equivalent to that of 
ordinary Yang-Mills theory at leading order in the large N expansion. At fixed-A, small L 
and a range of f3, the deformed theory maps to a two-dimensional theory with electric and 
magnetic (order and disorder) perturbations, analogs of which appear in planar spin-systems 
and statistical physics. We show that in this regime the deconfinement transition is driven by 
the competition between electric and magnetic perturbations in this two-dimensional theory. 
This appears to support the scenario proposed by Liao and Shuryak [1] regarding the magnetic 
component of the quark-gluon plasma at RHIC. 
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1. Introduction 

Through numerical simulations on the lattice [2, 3] and the experimental program at the 
Relativistic Heavy Ion Collider (RHIC) [4-7], we know that QCD has a high temperature 
deconfined quark-gluon plasma phase at temperatures above T c = /3" 1 170 MeV, where T c 
is parametrically of the order of the strong scale of the theory. Through lattice simulations, it 
is also known that the pure gauge sector of QCD, Yang- Mills theory, has a low temperature 
confined phase and a high temperature deconfined phase [8] . While symmetry and universality 
arguments are useful [9], to date, there is no direct continuum field theory technique to address 
most aspects of this transition due to its non-perturbative nature. 1 

Our goal in this paper is to make progress in understanding the microscopic mechanism 
driving the deconfinement transition of QCD and related theories, hopefully providing new 
insights into the structure of the quark-gluon plasma in the temperature region around T c . 

lr The deconfinement transition can also be studied by using strong coupling lattice models [10], however, 
it is not known how to extend this to phases continuously connected to the continuum. The existence of a 
deconfined phase (in continuum) can be established in perturbation theory [11]. 



- 1 - 



As reviewed in Ref. [12], this is one of the important problems concerning the physics of the 
nuclear collisions at RHIC. 

Recently, two new methods have been introduced for studying aspects of the deconfine- 
ment transition in a variety of gauge theories. The gauge/gravity correspondence, as realized 
by string theory, is a powerful tool for studying certain strongly coupled gauge theories [13-15]. 
The theories for which a semi-classical limit of string theory is useful usually differ from QCD 
in some way such as the existence of a non-decoupled KK-tower of states or an absence of 
asymptotic freedom. Regardless, this approach has the remarkable virtue of allowing one to 
do detailed calculations in a host of strongly coupled systems, many of which are plausibly 
in the same universality class as QCD or a QCD-like theory. For some recent applications to 
finite temperature properties, see [16-18]. A second approach for studying deconfinement was 
developed in [19-22], where one considers the large N limit of four dimensional SU(N) gauge 
theories compactified on S 3 x S 1 . For a small S 3 , the theory reduces to a matrix model, and 
there is a calculable deconfinement transition. In this second approach the large N limit is 
important for achieving the thermodynamic limit. Motivated by these two inspiring examples 
we pose the following questions: 

Can we find a calculable deconfinement transition in an asymptotically free and 
confining gauge theory by using field theory techniques? Is this even possible as 
the transition itself is non-perturbative? Can we give a simple physical picture of 
the mechanism behind the deconfinement transition? 

The small <S 3 x 5^ example provides an existence proof that finding calculable examples 
of deconfinement transitions in asymptotically free gauge theories is possible, at least at 
N = oo [20-22]. 

In this work we introduce field theory techniques through which the deconfinement tran- 
sition of four-dimensional Yang-Mills theory can be moved to a semi-classical domain where 
it becomes calculable using two-dimensional field theory. We achieve this by studying a 
double-trace deformation of Yang-Mills theory on R 2 x S[ x S^, which we refer to as the 
"f^-deformed Yang-Mills" theory or simply deformed Yang-Mills. Our deformation is similar 
to the one studied in the context of large N volume independence [23] . 

The O^-deformation has the effect that at large N, fixed-L, and arbitrary /3, the thermo- 
dynamics of the deformed theory is equivalent to that of ordinary Yang-Mills at leading order 
in the large N expansion. This thermal generalization of volume independence is depicted in 
Fig.l and described in Section 2. At finite N, thermal volume independence implies that the 
phase and thermal properties of the deformed theory in the interval: 

volume independence : L > (AiV) -1 (1-1) 

for a given j3 must coincide with the finite temperature properties of ordinary Yang-Mills 
theory up to 0(1/N 2 ) corrections. In the regime (1.1) the deformed theory remains incalcu- 
lable (without using lattice simulations) and the deconfinement transition cannot be studied 
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Figure 1: (Left) At N = oo, the thermodynamics of the deformed and undcformcd Yang-Mills 
theories are equivalent. (Right) Unlike pure Yang-Mills, the deformed theory has, at finite N, a semi- 
classical domain (defined by LA A < 1) where the confinement-deconfinement transition is analytically 
calculable. 



analytically. However at smaller L, a calculable, semi- classical regime opens up: 

semiclassical domain : L < (AN) -1 . (1-2) 

In this interval and at f3 = oo, we have analytic control over the three-dimensional long 
distance dynamics, and a weakly coupled (yet non-perturbative) semi-classical realization 
of confinement [23-25]. 2 For finite /3, the effective description of the thermal theory is a 
two dimensional system with electric and magnetic perturbations, in which /3 appears as a 
parameter. The confinement-deconfinement transition can be studied analytically within this 
2d model by varying f3, for arbitrary rank Yang-Mills theory, in contradistinction with [20]. 3 
The transition is plausibly smoothly connected to the deconfinement transition of the fully 
four-dimensional Yang-Mills theory. 4 

In this sense we find that the deconfinement transition of four dimensional Yang-Mills can 
be studied via a two dimensional field theory with electric and magnetic (order and disorder) 
perturbations. This two dimensional system has parallels to the statistical mechanical systems 

2 Such calculable regimes of QCD are, of course, not new. An analogous situation occurs in describ- 
ing hadronic matter at high density, where asymptotic densities provide a weak coupling (but again non- 
perturbative) calculable framework where expected features of hadronic matter at lower density is repro- 
duced [26,27]. Our deformation is, of course, more abstract, but morally similar. 

3 The thermodynamic limit is achieved without need for N = oo as our set-up is at infinite spatial volume. 

4 We emphasize that the existence of a semi-classical domain in the deformed Yang-Mills theory and the 
absence of such a domain in ordinary thermal Yang-Mills theory does not contradict thermal large- N equiva- 
lence, and volume independence. The semi-classical regime is LNA < 1, i.e., the scaling regime L ~ (7VA) -1 as 
N — > oo, while the domain of large- N equivalences is the regime L = O(N ) as N — > oo. As iV is increased, the 
semi-classical domain shrinks to a narrow sliver and the thermal large- N equivalence holds for any L — O(N ). 
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studied in [28]. In this framework the microscopic mechanism behind the deconfinement 
transition is a competition between electric and magnetic objects in a manifest and calculable 
way. At higher temperatures the electric objects are more relevant resulting in deconfinement, 
while at low temperatures magnetic objects are more relevant, resulting in confinement. 

The fact that the deconfinement transition is manifestly driven by a competition of electric 
and magnetic degrees of freedom we feel is the most interesting qualitative aspect of our work, 
and the one most prone to generalization to other four dimensional gauge theories, perhaps 
including QCD. In fact, the idea of deconfinement as a competition between electric versus 
magnetic objects has already been introduced into real world QCD by Liao and Shuryak as 
a possible way to explain some of the most interesting features of the quark gluon plasma 
at RHIC, in particular its relatively low viscosity to entropy ratio. Within the context of a 
simple toy model [1,30] (also see [31]) Liao and Shuryak argue that the viscosity to entropy 
ratio of a plasma of electric and magnetic excitations is minimized when the densities of 
magnetic and electric objects are comparable. 

The organization of this paper is the following. In Section 2 we elaborate on thermal 
large N volume independence and introduce the deformed Yang-Mills theory. In Section 3 
we discuss how the deformed Yang-Mills theory experiences weakly coupled confinement on 
R 3 x at small L. In Section 4 we describe the calculable deconfinement transition and the 
role of electric and magnetic objects. In Section 5 we conclude and discuss novel directions 
suggested by this work. 



2. Thermal large- N equivalence 

Consider ordinary four-dimensional Yang-Mills theory compactified on R 2 x T 2 , parameter- 
ized as: 

.M 4 (L,/3) = R 2 xSixSj, (2.1) 

where is a thermal circle of size j3 while is an ordinary circle of radius L. The action 
is: 



-yYM 



-/ 

J Ma 



^F^ + iO^rF^ 



(2.2) 



where F^ v = F^t 0, is non-Abelian field strength, F^ v = \^ vpa F pcr , g is 4d gauge coupling, 



and 9 is the theta angle. 5 For simplicity we henceforth set the 6*-angle to zero. 

This theory possesses a global (Zjv)/? x {%n)l center symmetry. This symmetry is the 
set of local SU(N) rotations periodic up to an element of the center group of SU(N): 



g(xi,x 2 ,x 3 + /3,x 4 ) = zp g(x 1 ,x 2 ,x 3 ,x 4: ), 
g(xi,x 2 ,x 3 ,X4 + L) = z L g(xi,X2,x 3 ,X4), 



z /3 — Z L — L - 



(2.3) 



moded out by the set of local gauge rotations (which are by definition single-valued on x 
Si). The order parameters for the (Zjv)/3 x {1^n)l center symmetry are the non-local Wilson 



3 We normalize the generators t a of the Lie algebra in the defining representation as: tr{t a t h ) — |5 a 
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lines: 



= exp[i / A^dx^], VLr = exp[i / A30IX3} (2-4) 

along the x circles, respectively. The center symmetry acts on the order parameters 
as: 

(Z N )f3 x (Z N ) L : trtyg -> trO^ 

: trSl L ^ z L trf) L (2.5) 

We define the "f^-deformed Yang-Mills" theory or simply "deformed Yang-Mills" as: 

[N/2] 



S dYM = S YM + AS, AS^ f V a n \tv(n n L ) 



(2.6) 



with sufficiently positive coefficients {a n } and [N/2\ denoting the integer part of iV/2. 6 
In the decompactifcation limit (3 — > 00, the deformed theory enjoys volume independence 

[23]: 

Large- N volume independence: Yang-Mills theory on R 4 is equivalent to the 
deformed Yang-Mills theory on R 3 x for any finite value of L, up to 1/N 2 
corrections, provided the (Z/v)l center symmetry is not spontaneously broken. 

Our deformed theory satisfies the condition that (Zjv)l remain unbroken by construction. 
Equivalent means that correlation functions of neutral sector observables - operators which 
are neutral under the (Ztv)l center symmetry - are the same in the two theories up to 
0(1/N 2 ) corrections. Volume independence does not apply to correlators containing non- 
neutral sector observables, the simplest example of which is trfii. Since many interesting 
physical observables are in the neutral sector, many interesting observables in pure YM theory 
can be extracted by studying the correlators in the deformed and volume reduced theory. 
The large- N volume independence theorem has an immediate generalization to general 

pe (0,00): 

Thermal large- N equivalence: Yang-Mills theory on R 3 x is equivalent to 
deformed Yang-Mills theory on R 3_fc x (S\) k x for any finite value of L and for 
a given /?, up to 1/N 2 corrections, provided that the [(Z^r)/^ center symmetry is 
not spontaneously broken. 

As before, in our deformed theory where k = 1 the condition that (Zjv)l remain unbroken is 
satisfied by construction. Also as before, the equivalence only applies to correlation functions 

6 Double-trace operators are also used in the works of Ogilvie et.al. to study phases with partial center 
symmetry breaking [42-44] . Pisarski and collaborators give a phenomenological effective theory description of 
deconfinement by using such deformations [40], and study various aspects in [41]. 
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of neutral sector observables. This important corollary to the volume independence theorem 
can be proven by a simple modification of the arguments of [23]. Instead of discussing the 
proof, we will state its main physical implications. 

Perhaps the most important part of the thermal equivalence is that the (Zjv)/3-center 
symmetry is a spectator symmetry. It is left intact during the projections and deformations 
which are used in the chain of equivalences that are used to prove the volume independence 
theorem. Thus, the expectation values and connected correlators of topologically non-trivial 
Polyakov loops (which are neutral under (Zjv)l, but charged under (Zjv)^) are also part of 
the neutral sector to which the thermal large- N equivalence applies. Thus these observables 
must agree in f^-deformed Yang-Mills and pure YM theory at any /3. This means, the 
thermodynamics of the two theories are part of their respective neutral sector dynamics. 

Thus, at leading order in N, the thermal Polyakov loops must agree, both in the confined 
and deconfined phases: 

<^tr^) dYM (L) = <"ti%) YM = J °' COn6ned (2.7) 

K N 131 V ; K N 131 \z p , /3</3 c deconfined 1 ; 

where zp is some iV-th root of unity. Furthermore the deconfinement temperature and the 
latent heat associated with the phase transition must also agree: 



ft™ = p™ 



1 + <^: 



n dYM _ n YM 
i Hi —Hi 



1 + 



(2.i 



These agreements hold in the strongly coupled LNA » 1 domain where volume independence 
applies. As a consequence of this volume independence, these quantities are independent of 
L at leading order in N J 

2.1 Why bother? 

A common criticism of large- N volume independence and in particular of deformation equiv- 
alences is that it maps a strongly coupled gauge theory to another strongly coupled gauge 
theory, neither of which is analytically calculable. So, why bother? 

It is true that in the strict N — > oo limit, neither of the equivalent pairs seem to be any 
easier. However, as we shall explain in the next two sections, at finite- N, the same deformation 
serves to engineer a semi-classical domain at small L. This semi- classical domain is continu- 
ously connected to the strongly coupled regime of the undeformed theory. Such a step usually 
cannot be achieved within the undeformed theory itself. What the deformation achieves is a 
generalization of Yang-Mills theory that depends smoothly on an extra-parameter. 8 



7 The matching of the deconfinement temperature is numerically tested in lattice regularized theory by 
simulating Yang-Mills theory with heavy adjoint fermions, a theory which emulates deformed Yang-Mills, on 
a [(1) 3 ] L x 2p [47], and agrees with large-scale lattice studies [8] . 

8 Small- volume or reduced formulation are also useful numerically. For example, QCD with adjoint fermions 
satisfies volume independence if one uses periodic boundary conditions for all fields [39]. Simulations and 
studies of the reduced QCD and related models appears in recent works [45-49] . 
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In this way one obtains a new weakly coupled regime of locally four dimensional gauge 
theories. Based on our experience with our best understood examples of non-perturbative 
quantum field theory, it is often useful to understand the various weakly coupled domains 
before attempting to understand the theory at strong coupling. Furthermore, the extension of 
calculations to the border of their validity can sometimes yield interesting information about 
the physics in the incalculable coupled domain. 

3. Weak coupling confinement 

It is easily seen that in the regime LNA < 1, the resulting effective long distance theory is 
a three-dimensional theory which enjoys weakly coupled confinement for a wide range of (3. 
We consider first the limit (3 — > oo which has been studied in [23]. We briefly review their 
derivations to establish the context and notations for the next section. 

In the zero temperature, weakly coupled domain, the deformed theory has a unique 
center-symmetric minimum for the Wilson line, Sl^. The fourth component of gauge field 
behaves as a compact adjoint Higgs field, and the theory reduces to a 3d Yang-Mills-Higgs 
system. The vacuum expectation value of the Wilson line is 

n L = n Diag (l, e 2 *" N , e^ N , . . . , e 2 ^" 1 )^) , (3.1) 

where r] = e m ^ N for even N and f] = 1 otherwise, up to conjugation by gauge rotations. This 
leads to Abelianization (or adjoint Higgsing): 

SU(N) -> U(l) N -\ (3.2) 

of the long distance dynamics. Since the fluctuations of eigenvalues are small due to the weak 
't Hooft coupling, the Abelianization holds quantum mechanically. 

Due to gauge symmetry breaking, the off-diagonal components of the gauge field acquire 
masses. The spectrum of the gauge fluctuations in perturbation theory is composed of levels 
each of which is TV-fold degenerate. The level spacing is j^. The masses and charges of the 
lightest VF-bosons are 

mw i = Jjji Qw t =gai, i = l,...,N. (3.3) 

Here, ai G A° ff — {a^} are the simple roots of the Lie algebra and ayv = — a i 1S the 

affine root (which is there due to compactness of the adjoint Higgs). A° ff is called the affine 
(extended) root system of the the associated Lie algebra, 

A aff = {"i>«2, • • • ,a N _i,a N } . (3.4) 

The roots £ A° ff obey: 9 

^ ■ aj = 5 it j - - , i,j = 1,...N . (3.5) 

9 We changed the normalization with respect to Ref. [23] for convenience, such that the simple roots nor- 
malize to unity. 
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Since the gauge symmetry is broken as SU(N) — > [/(l)^ -1 due to the compact Wilson 
line (3.1), there are N species of monopole-instantons. The topological and magnetic quantum 
numbers of these instantons are: 

and the negation of those for the anti-instantons. 10 

In three dimensions, Abelian duality relates a photon to a compact scalar. With cr J (x) 
the compact scalar dual to the photon -A^(x) of the j-th U(l) subgroup, the Abelian duality 
relation is: 

g 2 

F iiv = 7—7 <W d p cr . (3.7) 

To all orders in perturbation theory, (ignoring topological sectors), the long distance descrip- 
tion is free Maxwell theory in 3d, and is given by: 

The proliferation of these instantons leads to interaction terms in the Lagrangian of the 
compact scalar [32]. This is the generalization of Polyakov's mechanism to a locally four 
dimensional gauge theory [23]. The action for the low energy effective theory (the dual 
Lagrangian) in the small domain is 



6" 



dual 



N 

L (^) Z (V < T) 2 -C^cos(a l - < T) + .-. , 

i=l 



where £ = C e~ s ° = Am^ (g 2 N)~ 2 e - 87T2 /(9 2 N(m w )) j g ^ e monopole fugacity, and So is the 
instanton action. 

The action (3.9) is a non-renormalizable low energy effective theory valid at distances 
larger than ~ LN. Ellipsis stands for higher order terms in the semi-classical expansion 
as well as terms due to the omission of VF-bosons. 

The existence of a mass gap and linear confinement can easily be derived using the dual 
Lagrangian (3.9). The mass gap for the (N — 1) photon species is: 



/tt P \ t . ( LNA \ 5//6 

m„ = m a sin I — , m a = AK 



In 



/ LNA\ 



8/11 

\NJ' """ ""^V 2vr ) 



\ 2ir 



1,...,N-1 

(3.10) 



10 These monopoles are finite action topological configurations in the Euclidean formulation, and hence 
instantons [32] . N — 1 of them are ordinary 3d instantons, and the extra instanton, which has no counterpart 
in a microscopically 3d theory and which is pertinent to locally 4d nature of the theory, is sometimes called a 
twisted instanton. In a center symmetric background, all of these instantons carry equal action. Sometimes, 
they are also referred to as BPS and KK monopoles, or monopole-instantons. Some results about Bogomolny- 
Prasad-Sommerfield (self-duality) equations, and relation between these 3d and 4d instantons are reviewed in 
the Appendix. A. 
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where A is an 0(1) coefficient. The result of the semi-classical analysis is reliable in the 

LNA 

< 1 (3.11) 



2?r 

domain. 11 Obviously, the precise quantitative features of the mass spectrum and the string 
tensions have a non-trivial L dependence in the semi-classical domain. At strong coupling 
LNA 3> 2ir and at leading order in the large- N expansion, all the neutral sector observables 
must saturate to constants independent of L due to volume independence (see Section 2). 
Naturally, one expects the semi-classic description to match the strong coupling description 
around LNA ~ 1. 



4. Calculable deconfinement 

The deformed Yang-Mills theory on x R 3 exhibits weakly coupled confinement in the semi- 
classical domain (LNA < 1), where the theory experiences adjoint Higgsing. This "Higgsed" 
regime is analytically connected to the LNA 3> 1 regime and to the theory on R 4 in the sense 
that there exist no order parameters which can distinguish the two-regimes. In this section, 
we develop a formalism in the semi-classical domain which permits us to study the thermal 
phase transition. To do so, we consider a finite temperature compactification of the deformed 
theory on x R 3 which corresponds to the theory on Aii(L, (3) = x R 2 x at arbitrary 
P. 12 

First, let us momentarily ignore M^-bosons. (This assumption and its region of validity 
will be examined below.) At asymptotically low temperatures, ft 3> rn' 1 , the dynamics is 
that of the 3d dual theory (3.9). A more interesting regime is 



m 



} <C p < m^ 1 ~ m^e So/2 , (4.1) 



w 

where the size of the monopoles is much smaller than ft which in turn is much smaller than 
the inter-monopole separation. In this regime the potential induced by a monopole which is 
1/r in 3d is enhanced to log(r) at large distances which is the Coulomb potential of a charge 
in 2d. This can be seen by using the method of images from electrostatics. To incorporate 
this effect in field theory, it suffices to compactify the low energy effective theory (3.9) down to 
2d. In this domain, the theory reduces to a well-known two dimensional theory of "vortices" , 
which are the dimensional reduction of 3d instantons. The action is: 



^idual 



R 2 



N R 2 

L2 (V<t) 2 - Cm E cos («< • CT ) + • • •] ' a = Z (£) ( 4 - 2 ) 

i=l 



ii 



The power of logarithm, given in given in Eq.(3.37) of [23] as 9/11, is a minor error. More importantly, 
the small parameter in the the discussion of Ref. [23] is actually ), not LNA. Although the former is 

manifest in the formulae, the factor of 2tv was not explicitly written. It is actually useful to restore it. 

12 The analysis of the thermodynamics of the deformed Yang-Mills theory is analogous to the one of the 
finite temperature 3d Georgi-Glashow model [33-36] however, our set-up differs from it in the sense of being 
locally four dimensional. 
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where Cm = PC- For N = 2, this is the sine-Gordon model in d = 2 dimensions, and it is its 
generalization for N > 3. Whether a mass gap for the cr field is generated or not is tied with 
the question of the relevance of the e* ai <T operator. The conformal dimension of the operator 
e iai-a- a b ou t the free scalar fixed point is: 

Are**-! = °L = — = — • (4 3) 

1 1 4vra 4vra /3<? 2 " 1 j 

Note that for all oi; L 6 A° ff , the conformal dimensions are identical, because the algebra is 
simply-laced, af = 1, i = 1, . . . , iV. 13 

The perturbation of the free theory by the vortices is relevant if the conformal dimension 
A is less than two, irrelevant for A greater than two, and marginal otherwise. The quantum 
theory of (4.2) undergoes a phase transition at A = 2, 

An = ^ (4-4) 
T 

where subscript m stands for magnetic, between a phase of finite correlation length at low 
temperatures A < 2 (J3 > (3 m ), and a phase of infinite correlation length at high temperatures 
A > 2 (/3 < /3 m ). This is the well-known Berezinsky-Kosterlitz-Thouless (BKT) transition 
[37,38], albeit with an inverted temperature. In other words, the high temperature phase 
is populated by neutral magnetic vortex-anti-vortex pairs and these pairs dissociate at low 
temperature, opposite to the conventional BKT transition. This means, in the low and zero 
temperature phase, the mass gap is induced by the magnetic defects in the f^-deformed 
Yang-Mills theory. 

However, the effect described above is not the whole picture - the gapless phase is an 
artifact associated with the omission of electrically charged W-bosons, as noted in the context 
of the 3d Georgi-Glashow model in [34,35]. The W-bosons are not important in the long 
distance regime of the gauge theory on R 3 x because they are finite energy (mass) particles, 
as opposed to 3d instantons which are finite action defects. However, when the space is further 
compactified to R 2 x x S|, W-bosons traveling around the thermal circle have finite 
action, equal to f3mw- Their Boltzman weight is e~^ mw and has an interpretation as a W- 
boson fugacity. The W-bosons are a small perturbation (with respect to topological defects) 
when e~P mw <C e~ So or /3 3> Clearly, the scale at which monopoles become irrelevant 
is outside this regime. 



13 The ellipsis in (4.2) stand for perturbations sub-leading in the semi-classical expansion. Here, there are 
some subtle issues. Even at order k in the expansion, there is a sub-class of operators which has the same 

scaling dimension as the leading term, for example, A[e~ fe ? r " e l(Qi+ai + 1+ -' ai + fe) ' <T ] = A[e~ s 2 " e lai '"\ due to 
a Lie algebra identity, (at + aj+i + . . .a^+fc) 2 = 1 for k ^ N. In the effective theory, these and a plethora 
of others are there, generated and relevant in the sense of Wilsonian renormalization group. Although the 
scaling dimensions of these operators are identical to the ones that appeared in our effective Lagrangian (4.2), 

in the weak coupling domain, their prefactors are suppressed by extra-powers of e s 2jv . Hence, they remain 
as small perturbations at distances where the leading magnetic perturbation becomes strong. Thus, the effect 
of sub-leading terms are negligible there. Close to the boundary of semi-classical window, these operators may 
and will become important, as well as possibly near the critical temperature for the deconfinement transition. 
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p, 



A [cos 0^- o ] < 2 



m \ oo (3 



!p 

A [cos a. a ] < 2 1 * 

Figure 2: For /? > /3 m , e 4Qi <T is relevant, and for /3 < (3 e , e %ai ' a is relevant. In /3 m < (3 < (3 e interval, 
both perturbations are relevant. 



If one ignores the topological sectors of gauge theory, which is justified if e ^ mw 3> e 



(P ^ "g 2 ")' P ron fe ra tion of the two-dimensional gas of PV~-bosons generates an effective 
theory 

, ~ n i i r 

5 P.t. = / a 2 -( w ^ cos( ai • ?) + ... , a = — ^ = — 2 (4.5) 

where Ciy ~ ^e~^ mw and 4^ *dor = is the dual of cr field in 2d. The conformal dimension 
of the Ty-boson operator is 

• ~ a 2 Bci 2 
A[e ia ^] = -p- = 4vra = f^- , (4.6) 
47ra AttL 

The ellipsis in (4.5) stands for electric perturbations sub-leading in e~^ mw expansion, and 
the analog of the discussion in footnote (13) applies. This theory has a BKT transition at 
A[e ia ^\ = 2 or: 

Pe = ^ , (4.7) 

r 

where subscript e stands for electric. It has a gapped phase at high temperatures (3 < f3 e 
induced by free electrically charged excitations and a gapless phase at low temperatures where 
electrically charged excitations form neutral molecules. This makes sense because in the 
absence of topological defects, the large-/? theory is the compactification of the free Maxwell 
theory (3.8) which is related to the gapless phase of (4.5) via an a o 1/a or T-duality. 

The magnetic monopoles are a small perturbation (with respect to PF-bosons) when 
e — So ^ e -fim w or p ^ 4ti\L^ Clearly, the scale at which iy-bosons become irrelevant is 
outside this regime. This implies that neither electric nor magnetic BKT is actually there 
while (4.2) and (4.5) are valid descriptions. 14 



14 If, in Fig. 2, /3 m were larger than f3 e within the regions of validity of (4.2) and (4.5), this would have 
implied the presence of two genuine BKT transitions with an intermediate gapless phase. Indeed, in the planar 
Heisenberg model of ferromagnetism with a symmetry breaking perturbation (which reduces the symmetry of 
the theory to Zjv), such a phenomena takes place for all N > 4 [29]. 
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At arbitrary (3, and in particular, in a domain where both electric and magnetic pertur- 
bations are relevant, we should instead consider a Lagrangian of the form: 

. N N 

S= -(Vcr) 2 - Cm ^2 cos(ai- cr) - (w ^ cos(ai- a) + . . . , (4.8) 
•^ R2 i=i i=i 

where in the path integral we have to impose the duality relation * dtr = da as a con- 
straint. The electric-magnetic Coulomb gas representation associated with the field theory 
has the form V\ n t = Ve-e + V m - m + V e - m where V e - e (V m - m ) is the mutual logarithmic 
Coulomb interaction of electrically (magnetically) charged excitations and V e -m is the inter- 
action between electrically and magnetically charged excitations. For a detailed description 
and references to earlier related works, we recommend the reader Ref. [36]. 

4.1 Electric-magnetic competition and relation to Polyakov order parameter 

In the deformed Yang-Mills theory, the confinement-deconfinement transition is explicitly 
realized as a competition between electric and magnetic perturbations. This is a calculable 
realization of the scenario proposed in Ref. [1]. There are three regimes as a function of /3, 
as shown in Fig. 2. For /3 > f3 e , the e tai ' (T are relevant while the e iCti ' <T are irrelevant. In 
this phase, magnetic defects are free and dominate the long-distance dynamics, while the 
electrically charged particles are confined. For (3 < f3 m , the situation is reverted: the e iai ' a 
are irrelevant while the e iai ' a are relevant, which means that electric charges are free while 
magnetic defects are confined. In the interval j3 m < /3 < (3 e , both e iai ' a and e lCii ' a are relevant 
- we discuss this domain in more detail in Sec. 4. 2. 

In the small LNA domain, since the IR theory Abelianizes, the fundamental Polyakov 
loop may be identified with a " fundamental Quark" -operator. We define the following map- 
ping 

t 1 N 

—fig *~> — Ve^ (4.9) 

i=l 

where Vi, i = 1,...,N are the weights associated with the electric charges of the quarks in 
the fundamental representation. 15 If the external charge sourcing the Polyakov loop is in 



15 Conventions: The iV weights z/j are iV — 1 dimensional vectors forming an (N — l)-simplex. They satisfy 



Vi-Vi = \ tfy-^J, i,j = l,...,N. (4.10) 



(TV — l)-simplex is the figure associated with the defining representation of the algebra. At this stage, it is also 
useful to define the fundamental weights fik, 



Wb=$>j, k = l,...,N-l (4.11) 



Fundamental weights form the weight lattice A m , and the simple roots form the dual root lattice A r . A r is a 
sub-lattice of A w and the quotient is isomorphic to A w /A r — Zjv. The generators of the (A TO ,A r ) obey 

o i - W = i^ = i« iJ t,i = l,...,JV-l, (4.12) 
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some other representation, the mapping generalizes straightforwardly. 16 

The Lagrangian (4.5) apart from the obvious periodicity identification a ~ <x + 4iraj, 
is also invariant under a discrete 7Ln which we identify with the ordinary center symmetry, 
(Zjv)«. A shift in the weight lattice A„, acts as 

(Zjv)/3 : <x <r - 4-Kfi k 

In reaching the second step, we used the identities given in footnote. 15. Let us now calculate 
the the expectation value of the Polyakov loop. 

4.1.1 Low temperature 

In the (3 (3 e domain, we can safely use the Lagrangian (4.2) to describe the dynamics. The 
electric perturbations are highly suppressed and also irrelevant in the renormalization group 
sense. The insertion of an electric charge into the medium may be viewed as a vortex in the 
<r field theory. The vorticity is the electric charge associated with the probe, i.e., 

9. = J_ f da = u j (4.15) 

where C is a closed curve encircling the test charge. Thus, we have 

<t(<9) ~ UiO, V<r ~ -.eg (4.16) 

r 

To evaluate the action of a vortex in the free theory, we regularize the R 2 space to a disk 
D2(R) with radius R. (We also need a short distance cut-off. The finite size of vortex core 
serves this goal, but this short-distance divergence is unimportant for what follows.) It is: 

S(R)= [ d 2 x(Va) 2 ~[ rdr\~ log R (4.17) 
Jd 2 (r) J r 2 

This is, indeed, the Coulomb potential of a test charge in 2d and it clearly diverges as R — > oo. 
When we take into account the potential (4.2), we observe that the action grows quadratically: 
Sint.(R) ~ Id 2 (R) d 2 xcos(ai ■ a) ~ R 2 . However, this is an overestimation due to the form 
of the ansatz (4.16). The minimization of action in the space of possible <?{0) with the given 
vorticity generates a linearly rising action as R is increased. This is a configuration where cr 
is constant everywhere, but exhibits a jump along a cut. The punch-line is, in the confined 
phase, we have 



1 N 

i \ ^iui-cr\ Jim — - \ r I g — S[in the presence of i dcr=Vj}— 5[in its absence ] 

i- 



(4.18) 



the reciprocity relation. 

16 The generalizations of this mapping to anti-symmetric, symmetric and adjoint representations are: 

N N N 

trf^AS ~* J2 e i{Vl+V > y * trf^s — J2 e i{Vi+V i >B trfi^Adj — J2 e i{vt ~" s) '' . (4.13) 
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4.1.2 High temperature 

In the (3 <C (3 m domain, we can use the Lagrangian (4.5) reliably. The magnetic excitations 
are suppressed and irrelevant. Here we wish to calculate the expectation value of the e lUi ' a 
operator in a description where <x is the local field describing Lagrangian. The periodicity 
identification of the <x field is cr ~ a + Airotj. The potential V(<t) = —Cw ^2iLi cos(aj • <x) 
is also invariant under the (Zjv)/? center symmetry (4.14), and has N isolated minima within 
the unit-cell of the root lattice. This means that the theory has N thermal equilibrium states 
in this phase. The expectation value of Polyakov loop is: 

1 N 

(— J2 eiVi ' a ) = e+i ^ k > k = 0,l,...,N- 1 . (4.19) 

i=i 

These iV- minima can be rotated into each other by the action of spontaneously broken (Z^r)^ 
symmetry. 

This is precisely the picture that we believe should hold in Yang-Mills theory. In deformed 
Yang-Mills, we analytically demonstrated the existence of the two phases. 



Remark: The main novelty of this description is following: The effective dual Lagrangians, 
(3.9), (4.2) and (4.5), are already long-distance descriptions. The non-perturbative phenom- 
ena, such as a mass gap, linear confinement in the confined phase, and the existence of a 
deconfinement transition, are already in the tree-level description of the dual theory. This 
is a main difference between studies of deconfinement to date and our description. In our 
dual formulation, the long-ranged correlations are already built into the dual Lagrangians 
and correlation functions can be easily evaluated via these actions. This progress is possible 
because toroidal compactification with deformation introduces a new parameter, LNA, in 
the theory. When this parameter is taken large, we face the conventional problems of strong 
gauge dynamics. 

4.2 Estimate for phase transition scale 

In the semi-classical domain, the theory has at least two phases, (3 < f3 m where electric charges 
are free and magnetic charges are confined, and /3 > f3 e where magnetic charges are free and 
electric charges are confined. The phase transition must occur at some: 

P c e[Pm,Pe)= ^1 • (4.20) 

L 9 9 z i 

In this domain, we do not have a good tool to find the value of the transition temperature, 
as both perturbations are relevant. 

Despite the fact that we can demonstrate the existence of two phases (confined and 
deconfined) in a semi-classical approximation, the transition itself takes place in a regime 
(4.20) where the theory again becomes strongly coupled! 

We conjecture that the transition should occur when both electric and magnetic per- 
turbations simultaneously become order one following [36]. To argue this, note that if one 
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perturbation is order one while the other is small, then the system is gapped either due to 
electric excitations or magnetic excitations, which is to say the system is in one of the two 
phases. In such a case, the smaller effect may be treated within non-degenerate perturbation 
theory, and should not alter the behavior of the theory drastically. The two types of perturba- 
tions become comparable when the densities (i.e., fugacities) of electrically and magnetically 
charged quasi-particles become comparable. This is also argued to be the case in the scenario 

8tt 2 

of Ref. [1] within the context of QCD. Indeed, for SU(N), e ~^ mw = e~^Trk ~ e ~ s ° = e~~^ 
at: 

4vrL 4irLN 

Pc = — = , (4.21) 



9 2 A 

which is actually the midpoint of the (4.20). 

It is also instructive to study the conformal dimensions of the electric and magnetic 
perturbations around the deconfinement temperature. We find: 



A e = Afe 1 ^" 7 ] 



A R A r 



4vrL 

{Qw.Qku) 2 

(4vr) 2 



AttL' 



A, 



A[e 



4vr/3 g*/3 



1 



no sum over l 



(4.22) 



where the reciprocity of the dimensions of electric and magnetic perturbations is a consequence 
of the Dirac quantization condition 

— a<j = 47r (6ij — 2<5jj + i — \8i,j-i 



Qw, ■ QMi = 9&i 



9 



(4.23) 



At the critical point, the dimensions of both perturbations are equal to one, 



A[e* 



/3c 



A[e l 



1 . 



(4.24) 



In this sense, the theory as a function of /3 has four interesting domains and plausibly a single 



A m >2 


A m< 2 


A m <\ 


A m < 111 


A<><l/2 


A e <l 


\ A e <2 


A e > 2 oo p 



P. 



m 



p. 



p. 



Figure 3: A more refined version of Fig. 2. The putative phase transition is expected to occur at (3 C 
where A e = A m = 1. For j3 c < j3 < j3 e , magnetic operators dominate, and /3 m < j3 < j3 c , electric 
operators dominate. 



phase transition at /3 C , as depicted in Fig. 3. The theory exhibits confinement for f3 > (3 C , 
which corresponds to A m < 1, A e > 1 and deconfinement for j3 < (3 C which corresponds to 
A m > 1, A e < 1. 

This provides a more refined version of the domains of the thermal gauge theory relative 
to the Polyakov order parameter. We will speculate on the possible significance of f3 m in the 
conclusions. 
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4.3 Extrapolation to larger L or larger N 



The appearance of N in (4.21) is rather crucial, because the region of validity of the semi- 
classical analysis is LNA < 1, not LA < 1 (otherwise this would clash with large- TV" volume 
independence). At the boundary of the semi-classical domain, the transition temperature 
approaches /3 C ~ A" 1 , the expected result based on numerical simulations and dimensional 
analysis. For LNA 3> 1, by volume independence, this value must be saturated, up to 1/N 2 
corrections, hence the plateau shown in Fig. 4. The critical temperature is, T c = {3 ( 



-i. 



In 



' LNA^ 



LNA < 1 
LNA > 1 



(4.25) 

In order to reach to the volume independence domain, we do not necessarily need to increase 



< 
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odYM 

"ml.. 





1 2 
LNA 

Figure 4: Simplest possible phase diagram of SU(N) deformed Yang-Mills theory on R 2 x x S[. 
Above (below) the solid line, the theory is in the (Zjv)/3 unbroken (broken) confined (deconfmcd) 
phase. Between the solid and dashed line T c < T < 2T C , and at least in the semi-classical domain, 
although the theory is in deconfined phase, magnetic defects are still relevant. Below the dashed line, 
they are irrelevant. 



L. We can keep L-fixed while increasing N. The base space remains macro scopically two 
dimensional, but the dynamics (and thermodynamics) of the theory interpolates to that of 
Yang-Mills theory on R 3 x S^, In particular, the value of /3 C saturated in this regime must 
agree with ordinary Yang-Mills theory on R 3 x due to the finite temperature version of 
large- N equivalence. 

In Fig. 4, we plotted the simplest possible phase diagram of the theory. The semi-classical 
analysis is reliable in the LNA < 1 domain. We extrapolated the semi-classical result up to 
LNA/2ir = e" 1 ~ 0.367 where semi-classical function reaches to its local maximum. 7 

17 It is reasonable to ask up to what value of LNA one may expect that the semi-classical analysis will 
provide an accurate description. We guess (but do not have a solid argument) that this will be the case up to 
LNA ~ 2-7re~ 1 ~ 2.31. This question is in principle answerable by simulating deformed Yang-Mills theory and 
comparing it with our semi-classical results. 
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In the strongly coupled LNA 3> 1 domain, the transition temperature must be a constant 
due to large- N volume independence. Matching the transition temperature to the one of semi- 
classical analysis at the boundary of its region of validity, we obtained the phase diagram 
Fig. 4. It should be stated that in this phase diagram, the conjectural region is the vicinity 
of the matching point. Given T c at a value few times larger than the matching point, its 
L-independence at leading order in iV is dictated by volume independence. 

We also quote the numerical value for deconfinement temperature at LNA/2ir = e~ 1 
which we call the boundary of semi-classical window. If we use as the strong scale that of 
QCD, Ajyj-g = 213MeV, this gives us an estimate for SU(N) pure gauge theory 

T c = 0.74A « 158MeV (4.26) 

which is in the same ball-park with the lattice results, quoted in the Introduction. 18 

A final remark is in order for the SU(2) theory. Substituting a\ = —02 = 1 in (4.2) 
and (4.5), we observe that the discussion reduces to the one given in Ref. [34,36] for the 3d 
Georgi-Glashow model up to a trivial rescaling of the fugacity. (For SU(N) with N > 3, 
this is no longer the case, the analog of the monopole, which is on the same footing 
with a\, . . . , ajv-i does not exist in a locally 3d theory.) Ref. [34,36] argue that critical point 
resides in the (3 m < f3 < f3 e interval. They exhibit, by using fermionization, that the spectrum 
has a massless particle at criticality, and is of Ising universality. This agrees with universality 
arguments [9] and the numerical lattice studies for the SU(2) pure YM theory on R 3 x S^, 
and we expect the deconfinement transition to remain second order as the radius of is 
increased. (See Fig. 4). For the SU(N) case with N > 3, we were not able to determine the 
order of transition with confidence. We leave this for future work. 



5. Conclusions 

We introduced new techniques which enable us to continue the deconfinement transition of 
pure Yang-Mills theory to a calculable semiclassical domain. This was achieved by exploiting 
the recent developments in large- N volume independence and semi-classical confinement in 
gauge theories on R 3 x S\ [23-25]. 

Our approach uses a toroidal compactification of gauge theory on R 2 x x SL where 
at long-distances, the theory reduces to two dimensional field theory. A striking feature of 
this approach is that the deconfinement transition is manifestly seen to be due to a compe- 
tition between magnetic and electric perturbations in the two-dimensional field theory. At 
high temperatures the electric objects dominate, resulting in a deconfined phase. At low 
temperatures magnetic objects dominate, resulting in a confined phase. The order parameter 

Fixing the strong scale of Yang-Mills theories with n/ fermions, one finds, for n f = O(N ), that T c (n f ) = 
T c dYM /(l _ 2nj_y Numerically, these are T c {n f = 1) w 168 MeV, T c (n f = 2) « 180 MeV, T c (n f = 1) « 
193 MeV. Since for n/ > 1, the center symmetry is no longer an exact symmetry, the phase transition is 
replaced by a rapid-crossover. 
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distinguishing the two phases is the Polyakov loop, which is calculable in our framework away 
from the transition temperature, in both phases. 

The picture of the deconfinement transition as due to a competition between electric 
and magnetic objects, is a pleasing one. Confinement due to monopole instantons in 3d and 
due to a magnetic Higgs mechanism in 4d has a long precedent. More recently, Liao and 
Shuryak suggested the competition of electric and magnetic objects as a new way to look at 
the phase diagram of QCD. In our calculable deformation of Yang-Mills theory, their scenario 
is explicitly realized, at least in the semi- classical domain. 

In the semi-classical regime, both electric and magnetic perturbations are relevant in the 



window. In this regime, the densities of both electric and magnetic components of the plasma 
are comparable, while for T > 2T C the contribution of magnetic objects to the plasma is 
negligible. Could this window extrapolate to R 3 x of Yang-Mills theory? Below, we 
assume this logical possibility and speculate on its consequences. 

In Ref. [1] a plasma of electric and magnetic charges was studied using a classical molec- 
ular dynamics simulation with a variable electric to magnetic density ratio. The measured 
shear viscosity and diffusion constant were found to be lowest when the densities of electric 
and magnetic components are equal, and increased otherwise. Comparable densities of elec- 
tric and magnetic components arise naturally in our field theoretic description in the window 
(5.1). Our analysis implies that at T > 2T C , the magnetic perturbations become irrelevant. In 
this domain, magnetic excitations are confined to neutral molecules. If the model of Ref. [1] 
is a reasonable description of the relevant features of the quark-gluon plasma, and if we ex- 
trapolate our semi-classical results to the strong coupling domain, then we expect that for 
temperatures T > 2T C (which will be probed at ALICE detector at the Large Hadron Collider 
at CERN) a rapid increase of shear viscosity and diffusion constant with respect to RHIC 
results. 

5.1 Open problems 

There are many interesting directions that arise from our construction. Here, we sort a few 
which are most pertinent: 

1) It would be interesting to study the effect of fermionic matter on deconfinement in the semi- 
classical domain. In the presence of fermionic matter, the index theorem on R 3 x [50,51] 
implies that the mechanism of confinement is no longer necessarily due to simple monopoles, 
but rather due to magnetic bions and other non-self dual topological defects. (A classification 
of confinement mechanisms in semi-classical domain is given in [24].) 19 

19 The index theorem of Refs. [50, 51] is a variant of the well-known APS index for Dirac operator on 
R 3 x Si, a manifold with boundary. The index of a 4d instanton is a sum of indices of iV- types of 3d 
monopole instantons, Xid,inst. = X/i=i 3-M a . , each of which carry fractional topological charge 1/iV in the center 



2 



< T < 2T C 



equivalently 



Pm<p<P, 
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2) It would be useful to understand the order of phase transition in the NLA < 1 regime, and 
if possible, in the LNA ~ 1 domain for SU(N) with N > 3, both numerically and analytically. 

3) It would be interesting to generalize to orthogonal, symplectic and exceptional gauge 
groups, and in particular, to the groups for which the cover group has a trivial center sym- 
metry, such as G 2 - 

4) We have predictions for the L dependence of the critical temperature (3 C (L) (4.25), and mass 
gap (3.10) [23] in the semi-classical domain, and for their L-independence in the LNA ^> 1 
domain. It would be very interesting to test both regimes in lattice gauge theory and see up 
to what value of LNA the semi-classical description is in good agreement with lattice results, 
and at what value of LNA volume independence sets in. 

5) Supersymmetric gauge theories compactified on R 3 x are not expected to have any 
phase transition as a function of L [54] . A sub-class of supersymmetric theories such as pure 
M = 1 SYM, with gauge boson and adjoint fermion (A^,X a ), 20 also possess a semi-classical 
window in the LNA < 1 domain where confinement can be shown analytically. It would be 
useful to understand how deconfinement sets in when one consider this class of theories on 
R 2 x x SL with periodic boundary conditions for bosons and mixed 

X(x 1 ,x 2 ,x 3 + (3, x±) = - X(x 1 ,x 2 ,x 3 , x 4 ), 

X(xi, 22, X3, X4 + L) = + A(xi, X2, X3, 24). (5.2) 

boundary conditions for fermions. It may also be useful to understand how imposing periodic 
(supersymmetry preserving) boundary conditions in all directions avoids the phase transition. 

A. Yang-Mills in chiral basis and topological defects 

In this appendix, we remind the reader the topological defects pertinent to locally four dimen- 
sional gauge theories, in particular to R 3 x S[. It is useful to express the Yang-Mills action 
in a chiral basis which makes the role of self-duality manifest. We define r and the chiral field 
strengths F± (which furnishes (3, 1) © (1,3) representation of the Euclidean Lorentz group 
50(4) ~ SU(2) L x SU(2) R ) as 

Airi ~ i 

T — 2~ n ' -^*±/if — -^fj,u i -F^lu — F)j,v i . (A.l) 



symmetric background. A recent lattice simulation [52] calculates the index for some topological configurations 
and gives evidence for the existence of fractional topological charge objects. Outside the semi-classical window, 
the index theorem is still valid. However, since the topological defects are no longer dilute at large LNA or a 
large four-torus (which is suitable for lattice studies), it may be harder to probe fractional topological charge 
defects in this domain, see for example, [53]. Interestingly, the index obtained from lattice [52] and the one 
in [51] agrees. This is non-trivial and merits further study. 

20 Af = 1 mass deformation of J\f = 2 or Af — 4 SYM will work similarly. 
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The Yang-Mills action (2.2) can be rewritten as 



S YM = I i (r trF* - r tiF* ) (A.2) 



32vr v " +/» ""-/«/> 

The r — > ioo limit is the weak coupling limit. In the chiral basis, the instanton equation reads 

F_^ = or F tlv = F liV (A.3) 

For a 4d instanton, i*+ = &F 2 = 4F^F^, and topological charge is J tiF^F^ = 1. 
Its action and angle dependence appears as 

S = — / r trF? .„. = — / r tri^^L = 2vrir 

Thus, in the semi-classical expansion in 4d, the amplitude appears as 

e -S x = e -2nir = e -^T^ d ( A4 ) 

On small x R 3 , due to the center-symmetric Wilson line (3.1) associated with the 
boundary |x| — > oo, , there are more solutions to F^^ u + 0(g 2 ) = 0, Fij — heij^D^A^ + 
0(g 2 ) = 0, k = 1,2, 3, where 0(g 2 ) part, which is there due to deformation and one-loop 
potential, is omitted in weak coupling. The magnetic and topological charges of these A-types 
of monopole instantons are given in (3.6). In the semi-classical expansion, the amplitudes 
associated with these instantons are 

Notice that, the 4d instanton on R 3 x 5^ may be viewed as a composite of these N- types 
of 3d instantons associated with simple roots a% = 1, ... ,N — 1 and the twisted-instanton 
associated with affine root a at. These are sometimes referred to as "fractional instantons". 
The corresponding amplitudes obey 

N g7r a 

\{M ai = e~ 2 ™ = e -^~ ie . (A.6) 

8=1 

The semi-classical expansion on a center symmetric background is an expansion in ^ 2mr l N = 
e g 2 N n _ 'p^g 4^ instanton appears in this expansion at N order. In particular, at large- 
N, instantons are suppressed as e~ Sl = e~°( Nl ^ whereas the fractional instantons of center- 
symmetric background are e~ SMa i = e~°( N °\ hence they are part of large- N dynamics. 
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